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ABSTRACT 

MAD is the first multi-conjugated adaptive optics system at the VLT. We present H 
and Ks observations of the young massive cluster Trumpler 14 revealing the power of 
MCAO systems by providing a homogeneous Strehl ratio over a large field of view. 
Mean Strehl ratios of 6.0 and 5.9 per cent with maximum Strehl ratios of 9.8 and 12.6 
per cent in H and Ks, respectively, show significant improvement of the spatial PSF 
stability compared to single-conjugated adaptive optics systems. Photometry of our 
observations cover a dynamic range of ~ 10 mag including 2-3 times more sources 
than comparable seeing-limited observations. The colour-magnitude diagram reveals 
that the massive cluster originated in a recent starburst-like event 1 ± 0.5 Myr ago. 
We tentatively detect hints for an older population of 3 Myr suggesting that low 
intensity star formation might have been going on in the Hn region for a few Myr. We 
derive the luminosity function and mass function between 0.1 Mq and 3.2 Mq and 
identify a change of the power law slope of the mass function at m c ~ 0.53lg io Mq. 
The MF appears shallow with power law slopes of Ti = —0.50 ±0.11 above m c and 
T 2 = 0.63 ± 0.32 below m c . 

Key words: Instrumentation: adaptive optics - (Galaxy:) open clusters and associ- 
ations: individual: Trumpler 14 - Stars: pre-main sequence 



1 INTRODUCTION 

Adaptive optics (AO) systems have proven their outstand- 
ing ability in compensating atmospheric seeing over more 
than 20 years. With the first common user AO system, 
ADO NIS, mounted at the ESO 3.6m telescope on La 
Silla jRousset et al. 1990l L ground-based telescopes provide 
observations with previously unachieved angular resolu- 
tion. However, the performance of a single-conjugated AO 
(SCAO) system is restricted to small fields of view (FoV). 
With typical isoplanatic angles of 15 arcsec, the AO perfor- 
mance degrades quickly at distances of 20-30 arcsec from 
the guide star (GS). Thus, AO assisted observations are 
limited to maximum field sizes of ~ 1 arcmin with good 
AO corrections for a typical FoV of 30 arcsec. To provide 
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a spatially more stable performance, multi-conjugated AO 
(MCAO) systems are required. MCAO uses several GSs to 
correct the blurring due to Earth's atmosphere. In the frame- 
work of second generation instruments for the VLT and the 
E-ELT, the Multi-conjugate Adaptive optics Demonstrator 
(MAD), was developed in 2007 as the first MCAO system 
at the VLT (Marchetti et al. 2003, 2004, 2007). 

The proximity of the Carina Nebula and its cluster 
population provides an ideal testbed to assess the capa- 
bility of VLT-MAD. The combination of high spatial res- 
olution and the wide field of VLT-MAD allows to resolve 
the dense central regions of the clusters. NGC 3372, or Ca- 
rina Nebula, is a Galactic giant molecular cloud and a 
site of vigorous and ongoing star formation including sev- 
eral young massive clusters. It includes not only dozens of 
young O- and B-type stars but a lso evolved Wolf-Rayet 
stars fe.g. lMassev fc Johnson 1993l l. The youngest and most 
populous clusters of the region are the open clusters Trum- 
pler 14 and Trumpler 16 (hereafter Tr 14 and Trl6), lo- 
cated in the central part of the Carina Nebula. The high- 
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Table 1 . Strehl values over the lxl arcmin FoV 



Origin 


Filter 


max. Strehl 


min. Strehl 


mean Strehl 


Obs. MCAO 


H 


9.8 % 


3.3 % 


6.0 % 


Obs. MCAO 


K s 


12.6 % 


1.3 % 


5.9 % 


Sim. SCAO 


Hsim 


9.8 % 


< 0.01 % 


1.6 % 


Sim. SCAO 


K S ,sim 


12.6 % 


< 0.01 % 


2.0 % 



mass population with its strong UV radiation interacts with 
the surro unding material trigger ing subsequent star forma- 
tion (e.g. iRathborne et al. 2002T I. The stellar population of 
Tr 14 include s several O-type star s, notably the 02If* star 
HD93129Aa l|Walborn et al. 20021 ) which has been found to 
be the most massi ve star in Tr 14 wit h an estimated mass in 
excess of 100 Mq (|Nelan et al. 20041 '). Assuming a Kroupa- 
IMF, Tr 14 was sugge sted to be as massive as a few 10 3 
Mq (|Sana et al. 2010T ). which puts it close to the regime 
of the Galactic starb urst clusters like N GC 3603 YC or the 
Arches Cluster (e.g. iRochau et al. 20ld : IStolte et al. 20021 . 
resp.). The Carina Nebula is the closest region harboring 
such a massive young cluster, with recent studies estimat- 
ing the distance of Tr 14 to be between ~ 2 and 3 kpc (e.g . 
iTapia et al. 20031 ; ICarraro et al. 20041 ; lAscenso et al. 20071) . 
The lack of evolved stars makes Tr 14 a very young clus- 
ter. The ages of the high-mass conte nt of Tr 14 and Tr 16 
are estimated to be aroun d 1-2 Myr (|Vazquez et al. 19961 ) 
and 2-3 Myr fSmith 2006), respectively, while the formation 
of intermediate-mass stars probably started earlier. This is 
supported by the core-halo structure of the cl uster, with 
the halo being slightly older l|Ascenso et al. 20071) . However, 
the ages of the clusters as well as the distance estimates re- 
main controversial. One reason could be the anomalous ex- 
tinction law towards the clusters fe.g. ITapia et al. 20031 ). as 
this hampers precise age determinations, masking uncertain- 
ties/differences in the models for low- and high- mass stars 
in the main sequence and pre-main sequence, respectively 
l|Navlor 20091 ). 



2 VLT— MULTI-CONJUGATE ADAPTIVE 
OPTICS DEMONSTRATOR 

MAD has been developed as the first prototype MCAO sys- 
tem for the VLT. It aims at demonstrating the feasibility 
of MCAO systems in general and to investigate different re- 
construction techniques in particular. 

MAD is designed to correct for blurring due to atmo- 
spheric turbulence over 2 arcmin on the sky using Natu- 
ral Guide Stars (NGS). Two Deformable Mirrors are used 
for the MCAO correction. The first, optically conjugated at 
the telescope pupil, corrects for the ground layer turbulence, 
while the second, conjugated at 8.5 km elevation above the 
telescope, corrects for the field anisoplanatism. 

To examine the differences between the star- and layer- 
oriented reconstruction technique, two different wavefront 
sensors (WFS) are installed. The star-oriented MCAO cor- 
rection is supported by a Multi Shack-Hartmann WFS. The 
layer-oriented MCAO reconstruction uses a Layer Oriented 
Multi-Pyramid WFS. Both WFSs can sense simultaneously 
several NGS at visual wavelength but only one mode will 
be used at a time. Our observations were processed in star- 
oriented mode using 3 NGS. 

MAD is equipped with the IR camera CAMCAO 
(CAmera for MCAO) providing near-infrared broad- and 
narrow-band imaging, using a HAWAII2 2kx2k HgCdTe de- 
tector. With an image scale of 0.028 ar csec pixel" , the cam - 
era covers a 59 arcsec x 59 arcsec FoV |Amorim et al. 20061 ). 
A 2 arcmin x 2 arcmin FoV can be covered by moving the 
camera while the telescope is tracking and the AO loop re- 
mains closed. In contrast to SCAO, using multiple GSs al- 
lows to cover a larger AO corrected FoV with homogeneous 
Strehl ratios. 



We were able to obtain MCAO assisted observations 
with VLT-MAD to investigate the young massive star clus- 
ter Tr 14 (Pis: H. Sana, B. Rochau). In 2010, Sana et al. pre- 
sented their analysis of the data of the entire 2 arcmin x 2 ar- 
cmin FoV. Their scientific focus was on deriving the colour- 
magnitude diagram (CMD), structure of the cluster and to 
conduct a companion analysis. In contrast, our aim is the in- 
vestigation of the pre-main sequence population down to its 
low- and very-low mass members and the underlying mass 
function. Therefore, we concentrate on the deeper H- and 
-fCs-band observations of the central 1 arcmin 2 of the clus- 
ter. The results of our investigation are presented as follows: 
In Sect. 2 we introduce shortly the MCAO system MAD, 
followed by the description of the observations and the data 
reduction (Sect. 3). Sect. 4 includes a technical analysis of 
the MCAO performance. The scientific analysis is presented 
in Sect. 5 and our results are summarised in Sect. 6. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Observations 

Observations were carried out during the night of 2008, Jan- 
uary 12 by the VLT-MAD team. Deep H and Ks images 
were taken to map the innermost region of the young massive 
cluster Trl4 (a = 10 h 43 m 55 s , 8 = -59°33'03")- Short inte- 
gration times of 2 s were used to avoid severe saturation, and 
30 individual integrations are co-added to a single exposure 
of 1 min. Raw data reduction and image combination was 
performed using the ESO image processing software Eclipse 
l|Devillard 200 11 ). For sky and dark subtraction as well as for 
fiat field correction the calibration images taken during the 
science run are used. The sky and dark current subtracted, 
flat field corrected images are subsequently combined using 
the shift-and-add tool of Eclipse. 28 single exposures of 1 
min in each band are combined to a single frame of 28 min 
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Figure 1. Central region of Tr 14 as seen with VLT— MAD (north is up, east to the left). This HK$ colour composite shows the densest 
region of the cluster including the most massive cluster member HD93129Aa (a = 10 h 43 m 57.5 s , 6 = — 59°32'52.3" ) as the brightest star 
in the field. 
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Figure 2. Observed differences between the MAD and the NTT- 
Soft photometry after application of the constant ZP correction. 
ZP offsets are measured and plotted as a function of right ascen- 
sion (left panels) and declination (right panels) for the //-band 
(upper panels) and ifg-band (lower panels). Straight lines show 
the linear fit to the observed distribution revealing the spatial 
stability of the photometry. 



of total integration time. For comparison, the observations 
of the adjacent fields provided only 8 min of total exposure 
time (see also Table 2 of Sana et al. 2010). The final FoV of 
the central field covers an area of 68 arcsecx68 arcsec, and 
the HKs colour composite is shown in Fig.Q] Located close 
to the centre is HD93129Aa, the most massive member and 
brightest star in the field. 



3.2 Photometry 

To study the stellar content of Tr 14 , we perform PSF pho- 
tometry using IRAF/DAOPHOT (jStetson 199d ). Fitting 
15 stars in H- and Ks-b&nd, we find a Penny-PSF, com- 
prising a Gaussian kernel and Lorentzian wings, as the best- 
fitting PSF. It further turns out that a second order variable 
PSF provides a good match to the MAD data obtained with 
the star-oriented mode. Photometric calibration is based on 
NTT-Soft observations of Tr 14 published by Ascenso et al. 
(2007). The detectors and filter systems of MAD and NTT- 
Sofl are similar, leading to the absence of e.g. colour terms. 
Consequently, using the photometric catalogue of the MAD 
and NTT-Sofl observations is well suited to assess the good- 
ness of the PSF characterization of our MCAO observations. 

The comparison of seeing-limited NTT-Sofl photometry 
with our MAD photometry reveals no trend of photomet- 
ric zeropoint (ZP) with position on the MAD field (~1.1 
arcminxl.l arcmin) and shows a rms scatter of 0.11 mag 
in //-band and 0.12 mag in 7\s-band, respectively (Fig. [2]). 
Furthermore, Fig. [4] displays the ZP offset to be also con- 
stant with magnitude and colour. This agrees very well with 
the results of Sana et al. (2010) who applied PSF photom- 
etry to individual frames of 1 min instead of the combined 
image as used in our analysis. The spatially stable photom- 
etry emphasises that the stability of the PSF is comparable 
to that of seeing-limited observations, and exceeds the sta- 
bility of SCAO imaging, where anisoplanatic effects result 
in considerably larger variations in the PSF across the FoV. 

Thus, the combination of star-oriented MAD observa- 
tions with a second order variable PSF yields high preci- 
sion photometry across the entire FoV. The final photo- 
metric catalogue includes 1347 stars. Excluding stars with 
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Figure 3. _R"g-band completeness map of Trl4. Left Panel: i\s-band image of Tr 14, north is up and east to the left. Magnitudes at 
which 70 per cent completeness is achieved are superimposed as white contours (in steps of 1 mag). The image is centered (0/0) at 
a = 10 h 43 m 57.2 s , <5 = — 59°32'44.2" . Right Panel: i\s-band completeness map of the entire observed FoV. The grey shades correspond 
to the ii"s-band magnitude at which a completeness of 70 per cent is achieved. The strong local decrease of completeness close to bright 
stars without displaying large scale variations across the field shows that crowding effects are less important compared to the possible 
faintness of the stars. 
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Figure 4. Observed differences between the MAD and the NTT- 
Sofl photometry after application of the constant ZP correction 
as a function of magnitude (left panels) and colour (right panels) 
for the ff-band (upper panels) and i^g-band (lower panels). 

DAOPHOT photometric uncertainties above 0.1 mag (375 
sources), we used a catalogue of 972 stars for the further 
analysis comprising stars over a dynamic range of ~ 10 mag. 
Comparison to the 453 sources detected with NTT-SofI in 
the MAD field shows that VLT-MAD was able to detect 2-3 
times as many sources. 

3.3 Completeness 

We aim at deriving the luminosity function (LF) and, sub- 
sequently, the mass function (MF) of Tr 14. To derive the 
correct cluster LF and MF, accurate knowledge of the com- 



pleteness of the photometric catalogue is required. In dense 
stellar fields the detectability of a source depends on its 
flux and the local stellar density ('crowding'), and thus, on 
brightness and position of the source in the observed field. 
Crowding describes the effect of a decreasing detection prob- 
ability due to high stellar densities in the sense that a high 
brightness contrast between a point source and a possible 
bright neighbour limits its detectability. 

To determine the completeness, we apply the same tech- 
nique as described in Gennaro et al. (2011). We will briefly 
summarise the procedure here and refer to the Appendix 
of Gennaro et al. (2011) for further details. Artificial stars 
are inserted at random positions in the field, and the data 
analysis is re-run to estimate the fraction of recoveries. We 
added 50 artificial stars per run to our science image using 
the PSF derived with DAOPHOT. The rather low number 
of artificial stars has been chosen in order to not change the 
crowding characteristics of the observations. 100 individual 
runs are processed to add a total of 5,000 stars to our image 
in each photometric bin of 0.5 mag width. We achieve a typ- 
ical separation between simulated stars of d sim ~ 34 pixel. 

The final product of our procedure for each photometric 
band 'j' is the completeness as a function of position on the 
detector and magnitude of the star: 

^ , x a(x, y) . 

ex P( 7(*,y) ) + 1 

Here p is the magnitude of the corresponding star, a (< 1) 
is the normalization factor, /? the magnitude at which the 
completeness is a/2, and 7 describes how fast the complete- 
ness Cj drops to zero with decreasing brightness. With the 
derived values for (a, /3, 7), we can assign a completeness 
factor to each star considering its magnitude and position. 
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Figure 5. H- and _ft"g-band image of the core of Tr 14 (left and right panel, respectively). North is up and east is to the left. Overplotted 
contours show the level of equal Strehl ratios in H- and i\g-band. Apparent is the smooth decrease of Strehl values with the peak value 
centered on the GS included in the field in the A"s-band. The offset of the peak Strehl to the position of the GS in the i^-band is an 
artefact caused by the saturation of the star. Positions of the GSs are indicated by the grey diamonds and the encircled triangle by the 
straight grey lines. 



Fig. [3] shows the Tvs-band image of Tr 14 with the 70 
per cent completeness limits at different magnitudes super- 
imposed as contours (left panel). The right panel depicts 
the same completeness limits over the field as a function of 
magnitude. The limiting magnitude is significantly brighter 
in the proximity of the brightest stars. In less dense areas 
with fainter stars the completeness shows a smooth distri- 
bution over the field with a rather faint average limiting 
magnitude. In addition with the proximity of the cluster, 
this reveals that crowding is not a major source of incom- 
pleteness compared to the stellar faintness. Together with 
the largely increased number of detections when compared 
to the seeing-limited NTT-Sofl observations, this illustrates 
the improved performance provided by MAD. 



4 TECHNICAL ANALYSIS 

The observations have been carried out under observing con- 
ditions with a typical seeing between 0.7 and 1.6 arcsec. 
The MCAO correction provide data with an average full 
width at half maximum (FWHM) of 0.22 and 0.26 arcsec 
in H- and Ks-band, respectively. The relatively bad see- 
ing conditions, under which the observation were accom- 
plished, hampered a better correction. Considering the good 
geometry and brightness of the GSs together with good 
seeing conditions allow MAD to easily achieve AO correc- 
tions down to a FWHM of 0.15 arcsec, or even below (e.g. 
iMomanv et al. 2008; ICampbell et al. 2010l ). Seeing-limited 
NTT-Sofl observations of Tr 14 by Ascenso et al. (2007) 
had, for comparison, a FWHM around 0.7 arcsec. 

To investigate the performance of the MCAO system 



we measure the Strehl ratio of the combined image across 
the field. It is defined as the ratio of the measured peak in- 
tensity to the peak intensity of an ideal diffraction-limited 
image of a star of equal magnitude. Stellar properties are 
extracted from t he combined images usin g the Source Ex- 
tractor software (1 Bert in fc Arnouts 19961 ). The theoretical 
diffraction-limited PSF is created with the imgen task in 
the Eclipse package, using a primary and secondary mir- 
ror aperture of 8.2 and 1.1 m, respectively, a central wave- 
lengths of Ah = 1-65 fim and Xk s = 2.125 /im, with a 
width of AAh = 0.3 \im and AAk s = 0.35 /im. We derive 
the Strehl ratios of our H- and .Kg-band observations and 
the distributions of Strehl ratios in H- and ivs-band over 
the observed FoV are shown in Fig. [5] as white contours. As 
expected in the star-oriented mode of MAD, the Strehl ratio 
is maximised on the GS. Maximum Strehl values of 9.8 per 
cent (//-band) and 12.6 per cent (.Ks-band) are measured. 
The Strehl map, furthermore, reveals a shallow decrease of 
the Strehl ratio over the FoV. This is indicated by mean 
Strehl ratios of 6.0 per cent for //-band and 5.9 per cent for 
Ks-band. The results are summarised in TableQ] A steep de- 
crease of the AO performance, as observed in SCAO obser- 
vations, is alleviated due to the decreasing distance towards 
the other two GSs which are located outside our FoV. 

To quantify the difference between MCAO and SCAO, 
we compare our observations to simulated SCAO observa- 
tions. We create a grid representing our observed MAD FoV 
and assign values of Strehl ratios to each grid element with 
maximum Strehl ratios at the GS location according to our 
observations. The spatial variation of th e Strehl ratio is as - 
sumed to follow SR = S exp(-[J^] 5/3 ) l|Cresci et al. 20051 ) 
with So as the maximum Strehl ratio and an isoplanatic 
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Figure 6. Kg vs. H — Ks CMD of the central 0.75 pc of Tr 14. The cluster population can be identified as the well populated PMS in 
the CMD and goes down to (and beyond) the H-burning limit. Left Panel: The CMD is well represented by a Siess isochrone of 1 Myr 
in age (red solid line); red dashed lines show the uncertainties in the distance modulus. The MS isochrone is shown by the thin and the 
dotted black line. The box (dotted lines) displays the region that is shown in the right panel. Right Panel: The CMD with the 0.5 Myr, 1 
Myr and 3 Myr Siess isochrones overplotted (blue dashed, red solid and green dash-dotted line, respectively). Labelled masses are taken 
from the 1 Myr Siess isochrone. 



angle ©i so = 15 arcsec as a typical value for SCAO obser- 
vations. Results of the simulations are also summarised in 
Table[T] We found mean Strehl values of 1.6 and 2.0 per 
cent for the H- and Tfs-band, respectively. The larger mean 
Strehls in our MAD observations reveal the improved spatial 
stability of the PSF in our MCAO observations compared 
to SCAO imaging. 

Although SCAO observations provide higher peak 
Strehl ratios, the fastly decreasing performance with dis- 
tance to the GS makes a good characterization of the PSF 
at larger distances difficult. Especially in crowding-limited 
regions, such as dense stellar clusters, this hampers the cen- 
sus of the stellar population. The greater homogeneity of 
the PSF provided by MCAO observations allows, therefore, 
the coverage of larger AO corrected areas on the sky, crucial 
for photometric and astrometric studies of extended stellar 
populations. 



5 SCIENTIFIC ANALYSIS 
5.1 Colour— magnitude diagram 

Using the photometry of Sect . 13.21 we create a Ks vs. H — Ks 
CMD, shown in Fig.[6] The CMD allows the identification 
of the cluster population, primarily the cluster pre-main se- 
quence (PMS) population between 10 mag< Ks <16 mag 
and 0.2 mag< H — Ks <0.5 mag. The transition region 
from PMS to the main sequence (MS) is sparsely popu- 
lated but can be identified at 10.5 mag< Ks < 12.5 mag 
with 0.1 mag < H — Ks < 0.3 mag. Stars redder than the 
PMS (H — Ks >0.5 mag) suffer either from severe differen- 
tial extinction in the region, are background objects, highly 
reddened by the Carina Nebula or belong to the cluster as 
PMS stars surrounded by circumstellar material such as a 
disc causing the iCs-band excess. 



Without a comparison field, we cannot precisely quan- 
tify the contaminating fore- and background field star pop- 
ulation and, hence, subtract statistically the contaminat- 
ing field. However, the identification of different features of 
the cluster in the CMD (PMS, PMS-MS transition region) 
still allows the comparison with stellar evolutionary mod - 
els. For the PMS, we use Siess tracks (|Siess et al. 2 000). 
computed as in Da Rio et al. (2009), using the BT-Settl 
Grid of Allard et al. (2010). MS isochrones are from Marigo 
et al. (2008). Isochrones that resemble best the observed 
features are shown in the right panel of Fig. [6] superim- 
posed onto the CMD. The shape and width of the PMS- 
MS transition has proven to be a good a ge estimator (e.g. 
iBrandner et al. 20081 ; iRochau et al. 20101 ), but in the case 
of Trl4 the region is sparsely populated. Hence, to estimate 
the age of the cluster we use the colour difference between 
the MS and the PMS as well as the shape of the PMS. We 
identify the cluster PMS as best represented by the 1 Myr 
isochrone indicating that Tr 14 is a very young stellar pop- 
ulation. The derived distance modulus (DM = 11.8 ± 0.4 
mag) corresponds to a distance of ~ 2.3 ± 0.4 kpc. The 
uncertainty is estimated by finding the closest and far- 
thest distance which gives an acceptable approximation of 
the PMS of the cluster. The extinction was simultaneously 
derived to Ak k = 0.38 ± 0.03 mag, using Rv = 4.16 
l|Carraro et al. 2 004) and the relation of visual and near- 
infrared extinction according to Cardelli et al. (1989). The 
error on the average extinction is estimated by shifting the 
isochrone towards the blue and red boundary of the upper 
PMS, respectively. 

The distance of 2.3 kpc is comparable, though slightly 
smaller, to earlier derived distance estimates. For exam- 
ple, Carraro et al. (2004) locates Tr 14 at a distance of 2.5 
kpc, while Tapia et al. (2003) estimated a somewhat larger 
distance of 2.63 kpc. Comparing our extinction estimate, 
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Figure 7. Left Panel: Hess diagram based on our photometric catalogue. To distinguish the PMS stars without infrared excess from 
the contaminating field population and excess sources in the CMD we cut the Hess diagram into stripes of 0.5 mag in magnitude. Blue 
diamonds show the central magnitude of each stripe (y-axis) and the colour of the peak of the underlying stellar density (x-axis). The 
errors bars show the ltr-deviation from the fitted Gaussian function. PMS stars without significant infrared excess are identified by being 
inside the 1.5cr limit. Inlay Panel: An example of the density distribution for one magnitude bin (here: 14.6 mag< Kg < 15.1 mag). 
The black line shows the stellar density as a function of colour. A Gaussian function is fitted around the peak of the distribution and 
overplotted as the green line. Right Panel: The CMD as in Fig. [6] with the identified field contaminants and 'excess' sources highlighted 
(crosses and grey dots, respectively). The field contamination has been excluded while deriving the LF of Tr 14 while the 'excess' 
population is considered. Dotted lines mark the region of the excluded MS population. The dashed line is the 1 Myr MS isochrone by 
Marigo et al. (2008) and nicely displays the good identification of the field contaminants. 



which corresponds to an Av = 3.0 ± 0.4 mag following 
the applied extinction law, reveals a slightly higher aver- 
age extinction toward s Tr 14 compared t o earlier values, e.g. 
A v = 2.6 ± 0.2 mag ilTapia et al. 20031 ) or A v = 2.0 ± 0.2 
mag (|Carraro et al. 20041) . We benefit from a higher angular 
resolution, and the improved photometry might explain the 
slight differences observed in the different studies. 

The best fitting values, distance^ 2.3 ± 0.4 kpc and 
A Ks = 0.38 ± 0.03 mag, are derived using a 1 Myr PMS 
isochrone. This isochrone fits the PMS-MS colour width 
and shape of the PMS best. With a different combina- 
tion of distance and extinction (distance= 2.6 ± 0.4 kpc, 
Ak s = 0.36 ± 0.03 mag), a slightly younger isochrone (0.5 
Myr) represents the PMS-MS colour difference well, but the 
PMS shape is less well reproduced (Fig.|S}. For the 2 Myr 
isochrone the colour difference between MS and PMS is very 
narrow compared to the observations. We therefore adopt 
an uncertainty of 0.5 Myr on the age. However, very young 
stellar ages should be dealt with carefully as evolutionary 
models have large uncertainties for ages < 1 Myr, given the 
strong dependence of the predicted luminosities and tem- 
peratures on the in itial conditions of the computation (e.g. 
iBaraffe et al. 20021 ). 

Comparing the CMD with the Siess isochrones, we find 
the 1 Myr isochrone to best resemble the PMS. However, 
at 14 mag< Ks < f6 mag we identify a number of stars 
with colours bluer than the isochrone. Apart from observa- 
tional biases such as photometric uncertainties or a possible 
lower extinction of the identified sources compared to the 
derived average extinction, the observed sequence could al- 
ternatively be explained by a somewhat older population 



of PMS stars. By superposition, we observe that a 3 Myr 
isochrone nicely follows this blue sequence (see Fig. [6}. An 
older population would support earlier statements of con- 
tinuous star formation in the Carina Nebula, creating a 
surrounding older h alo population with ages up to 5 Myr 
l|Ascenso et al. 20071 ). 

5.2 Luminosity function 

5.2.1 Selection of different populations 

Since the observations have been carried out without ob- 
serving a comparison field, the isolation of the cluster popu- 
lation remained challenging, fn order to get a better idea of 
the cluster population we apply the following approach: We 
create a Hess diagram based on our H- and ifs-band obser- 
vations (left panel of Fig. [TJl which shows the stellar density 
in the colour-magnitude plane rather than the single stars. 

The Hess diagram is divided into 20 stripes, 0.5 mag 
wide, fn each stripe we examine the stellar density as a func- 
tion of colour. An example of such a distribution is shown 
in the inlay diagram in Fig. [7] Around each peak we find 
a rather symmetric distribution which is well represented 
by a Gaussian function (peak positions are highlighted as 
blue diamonds in Fig. |7J. The best- fitting Gaussian func- 
tion is shown by the green line superimposed to the stellar 
density distribution, ler deviations from the peak are plot- 
ted as horizontal error bars. The sparsely populated cluster 
MS yields unreliable fits. Hence, we excluded sources with 
K s < 12 mag and H-K s < 0.2 mag (dotted lines in Fig.0. 
This does not change our results, as we intend to focus in 
the present work on the PMS population. Thus, we aim to 
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Figure 8. _R",g-band LF of Tr 14. Displayed is the number of stars found in magnitude bins of 0.5 mag normalised to the observed area. 
The histogram shows the number of detected sources while the solid line depicts the same number but corrected for incompleteness. 
The dashed histogram is the completeness corrected LF including also 'excess' sources. The errors are the Poisson errors of each bin. 
Corresponding masses, highlighted by the arrows are taken from the 1 Myr Siess isochrone. The LF shows an increasing number of stars 
with decreasing magnitude to a turnover observed in the 0.25 — 0.54 Mq bin (PMS) and/or in the 0.18 — 0.25 Mq bin (PMS+'excess' 
sources). The dash-dotted lines mark the limit at which the completeness reaches 50 per cent and 70 per cent, respectively. 



isolate the PMS stars without infrared excess from lower 
mass MS stars belonging to the region but not to the young 
cluster and foreground sources, for simplicity in the follow- 
ing together called field contamination, as well as from the 
sources with .Ks-band excess. Comparison of the Hess dia- 
gram with the CMD reveals that a limiting deviation of 1.5a 
divides best the different populations. Sources being off the 
peak by more than 1.5a are considered as field contamina- 
tion (bluer than peak) or as 'excess' cluster sources (redder 
than the peak) and displayed by the crosses and grey dots 
in the right panel of Fig. [7] respectively. The limit has been 
chosen such that stars which can be clearly assigned to the 
contaminating field are separated from the PMS of the clus- 
ter. The PMS is nicely confined along a relatively narrow 
sequence and the chosen limit provides a good division of 
the PMS sources into a non-excess and an excess sequence. 



Although we cannot clearly identify the origin of the ob- 
served excess, these stars are considered as excess sources of 
the cluster. Firstly, the cloud of the Carina Nebula is effec- 
tively blocking contaminating background objects. Further- 
more, the effect of differential extinction should be lowered 
by the presence of the massive stars which removed the gas 
remaining from the star formation process such that highly 
reddened sources are unlikely to suffer from a large extinc- 
tion. We identified 79 sources with infrared excess and 274 
PMS stars without excess and, assuming all these stars to 
belong to the cluster, the fraction of excess sources would 
only amount to ~ 22.5%, which is lower than the excess 
fraction of 35% as reported by Ascenso et al. (2007). Con- 
straining the different regions of the CMD is in our case the 
best way to get a least biased view of the stellar population 
of Tr 14 and, subsequently, the best cluster LF. 



5.2.2 Ks-band Luminosity Function 

The LF is the number of stars as a function of magnitude. 
Stars are counted in magnitude bins of 0.5 mag, the counts 
are corrected for incompleteness and normalised to the con- 
sidered area. The final A"s-band LF is shown in Fig. [5] The 
histogram displays the LF based on stars that belong to 
the PMS stars without excess (black dots in Fig. [TJl and is 
not yet corrected for incompleteness. The completeness cor- 
rected LF is shown by the straight line with error bars that 
correspond to Poisson errors. The steep increase for stars 
fainter than Ks ~ 16.5 mag in the completeness corrected 
LF is due to the high incompleteness at the faint end. The 
dotted line further displays the LF of PMS stars including 
excess sources (black and grey dots in Fig. and is cor- 
rected for incompleteness. 

The LFs depict an increasing number of stars towards 
fainter magnitudes at the bright end until it turns over 
at K s = 14.75 mag (PMS-LF) or Ks = 15.25 mag 
('PMS+excess'-LF), respectively. This same turnover is also 
visible at Ks ~ 15 mag in the two Ks~ LFs derived by As- 
censo et al. (2007) based on NTT-Sofl and VLT-NaCo ob- 
servations, respectively. 

Similar AT-band LFs have been identified in clus- 
ters, observationally an d by simulating model LFs of 
synthetic clusters (e.g. IZinnecker fc McCaughrean 199ll ; 
lLada fc Lada 20031 '). Zinnecker & McCaughrean (1991) sim- 
ulated primarily if-band LFs for clusters of non-accreting 
PMS stars. They found time-dependent LFs exhibiting max- 
ima that shift with time. Comparison of our LF with its max- 
imum at Ks ~ 14.75 mag with the results of Zinnecker & 
McCaughrean (1991) reveal comparable functions for clus- 
ters between 0.7 and 1 Myr, consistent with the derived age 
for Tr 14. Such maxima are most likely caused by Deuterium- 
burning stars that lead to a peak in the LF. In a 1 Myr old 
cluster a PMS star of 0.3 Mq starts burning Deuterium, in 
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agreement with the masses that correspond to the observed 
peak in the LF. 

5.3 Mass function 

5.3.1 The Trumpler 14 Mass Function 

The initial mass function (IMF) represents the distribution 
of stellar masses for a given stellar system at the time of 
its formation. The form of the MF changes with time due 
to stellar, as well as dynamical evolution. As a result of 
such evolutionary effe cts, which affect entire stella r systems 
such as star clusters |Baumgardt fc Makino 20031 ). the MF 
is changed and, thus, the 'actual' MF will be called present- 
day mass function (PDMF). The PDMF is described, similar 
to the IMF, as a power law: 

d log N (log m) 

1 pd = r-j , (2J 

dlogm 

where iV(logm) is the PDMF. This slope is given by the lin- 
ear relation between the logarithmic mass intervals and the 
corresponding number of stars (in logarithmic scale). With 
this notation t he Salpeter-IMF has a slope of T = —1.35 
[jSalpeter 19551 ). Due to the youth of Tr 14, most of its mem- 
bers are populating the PMS which ends with the ignition of 
hydrogen in stellar cores at the MS turn-on. To obtain the 
mass of all PMS stars, we use the relation of the stellar lu- 
minosity and mass, the mass-luminosity relation (MLR), of 
the 1 Myr Siess isochrone. To convert apparent to absolute 
magnitudes we consider the distance of 2.3 ± 0.4 kpc and a 
constant foreground extinction of Ak s ~ 0.38 ± 0.03 mag 
as derived in Sect. 15. ll The mass of each star is estimated 
with the derived MLR and the corresponding magnitude. In 
case of the PMS stars without excess we have considered the 
brightness of the stars in Tfg-band, as this minimises the ef- 
fect of differential extinction. PMS stars with circumstellar 
material are, however, exhibiting an excess in the JCs-band. 
Therefore, we have used the i/-band instead to reduce the 
effect of the excess when deriving the mass. 

To obtain the PDMF of Tr 14, we use the method sug- 
gested by Mafz Apellaniz & Ubeda (2005). They discuss the 
biases that are i ntroduced by the b inning process to derive 
the MF (see also lMafz Apellaniz 20091 ). Assuming a constant 
bin size can lead to misleading results, due to the correla- 
tion between the number of stars per bin (higher for lower 
masses) and the assigned weights (from the Poisson statis- 
tics). To circumvent this problem they suggested a variable 
bin size to have an equal number of sources in each bin. They 
found, in this case, that the bias is almost independent on 
the assumed number of stars per bin, and remains low even 
for a single star in each bin. 

To minimise the effect of the binning and to maintain 
the statistical information, we have chosen to derive the 
PDMF in single star bins. Therefore, we sort the stars of 
our catalogue by mass. Each 'bin' consists of a single star 
and its bin width is defined to be the mean of the star and 
the preceding and following star, respectively. For the most 
and least massive stars in the sample, the upper and lower 
boundaries of the bins are chosen to be symmetric from the 
actual mass of the star. The 'number' (1 by definition) of 
stars is corrected for incompleteness associated with the star 
and divided by the width of the bin. 



The resulting PDMF is illustrated in Fig. [9] It displays 
two functions composed of i) only the PMS stars (index: 
PMSi) and of ii) PMS stars and the 'excess' sources (index: 
alL). In both cases, we identify a PDMF with a turnover 
at log(m) ~ —0.25. The jump at log(m) ~ —0.6 is due to a 
significant change in the slope of the MLR (~ 0.06 MQ/mag 
for m < 0.25 M and ~ 0.76 M Q /mag for m > 0.25 Mq). 
Therefore, and due to a probably increasing contribution of 
the Galactic field stars, the PDMF below log(m) < —0.6 is 
not considered in the following fitting process. Thus, we fit 
a two-component power law with a variable characteristic 
mass in the mass range 0.25 < M < 3.2Mq. 

i) PDMF of PMS stars without excess sources (black 
dots in Fig.^ty: 

We find a characteristic mass of m c = 0.52^o'i7 Mq 
and power law slopes of Fpd.pmsi = —0.46 ± 0.20 above ra c 
and rpD,PMS2 = 0.73 ± 0.59 below m c . 

ii) PDMF of PMS stars with excess sources (black and 
grey dots in Fig^ty: 

We derived power law slopes of FpD, a iii = —0.50 ± 0.11 
and TpD.aiw = 0.63 ± 0.32 above and below a turnover 
mass of m c = 0.53tg;io 

Mq, respectively. The slopes 
in both regimes agree very well within the fitting uncer- 
tainties. Consequently, no systematic bias is introduced 
due to the inclusion of excess sources and the deriva- 
tion of the PDMF slope is independent of our cluster 
member selection. Therefore, when addressing the PDMF 
of Tr 14, we refer in the following to the PDMF includ- 
ing the excess sources, i.e. to r P D, a iii = Ti = —0.50 ± 0.11 
and FpD,aii2 = ^2 — 0.63 ± 0.32 above and below m c = 
0.53+8:3 M . 

The turnover mass that is identified agrees very well 
with th e mass at whi ch the flattening of the Kroupa-IMF 
occurs (|Kroupa 200 ll ). Such flattening of the MF around 
~ 0.53 Mq is not an unique f eature of Tr 14, but observed 
in other clusters (e.g. Fig. 3 in lBastian et al. 20ld ). 

The slopes of the PDMF in the different mass regimes 
(above and below m c , respectively) are notably shallower 
than the corresponding part of the Kroupa-IMF (Salpeter 
slope r = —1.35). However, a source of uncertainty in the 
derivation of the MF slope are unresolved binaries. Depend- 
ing on the binary fraction and the intrinsic slope of the stel- 
lar MF, the observed slope may d eviate from the intrin- 
sic slope by u p to Ar ~ +0.5 (e.g. ISagar fc Richtler 199ll : 
iKroupa 200 ll ). Even with such a large correction our PDMF 
would appear natter than a Kroupa-IMF. 

Considering the rather flat slope of the PDMF of 
the central region of Tr 14 and the unaccounted effect 
of unresolved binaries, the derived PDMF is in good 
agreement with the central PDMFs observed in other 
young and massive clusters. In the very central region 
in NGC3603YC slopes of the MF are found between 
r ~ —0.31 and -0.91, depending on the stellar popula- 
tion (MS, PMS, P MS plus MS) and area considered in de- 
riving the MF fe.g.lSung fc Bessell 20041 : IStolte et al. 20061 : 
lHaravama et al. 20081 ). The different considered mass ranges 
do not hamper the comparison as above 0.5 M^ the 
MF is descr ibed by a single-power law (|Salpeter 19551 : 
IKroupa 200ll ). In the case of Westerlundl, Gennaro et al. 
(2011) found a significant flattening of the PDMF towards 
the centre of the cluster with slopes down to T ~ —0.7 al- 
though the very centre could not be considered due to severe 
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Figure 9. MF of Tr 14 derived following the suggested method of Mafz Apellaniz & Ubeda (2005). Each 'bin' contains a sin- 
gle star, which has been corrected for incompleteness and is divided by the width of each bin. The two MFs are composed of 
stars of only the PMS (left panel) and of the PMS plus the 'excess' sources (right panel). The best fits to the MFs are shown as 
straight red lines while dashed lines display the la interval. The 'high-mass' part is well represented by a power law with a slope 
of rpo pmsi = —0.46 ± 0.20 (PMS) and PpD alll = ~ 0.50 ± 0.11 (PMS+ 'excess' sources). The 'intermediate' part appears increas- 
ing with stellar mass (Ppd,pmS2 = 0.73 ± 0.59 (PMS) and rpD, a ll2 = 0.63 ± 0.32 (PMS+ 'excess' sources) with a turnover observed 
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saturation effects. In the central 0.75 pc, as in our study of 
Tr 14, Brandner et al. (2008) could identify a similarly flat 
PDMF with T ~ -0.5. 

However, it should be noted that we have to regard the 
effect of dynamical mass segregation in case of NGC 3603 YC 
and Westerlund 1 flattening the central PDMF because of 
the significantly higher density in these clusters. In contrast, 
no mass segregation was observed in the centre of Tr 14 in 
the considered mass range. Based on a minimum spanning 
tree analysis, mass se gregation was ide ntified for stars more 
massive than 10 Mq (|Sana et al. 2010l ). They have analyzed 
the same MCAO data of Tr 14 but including also the shal- 
lower MAD observations of the adjacent fields, such that the 
analyzed FoV encloses the central 2 arcmin. A comparable 
result has be en reported for the ONC with mass seg regation 
above 5 M l|Hillenbrand 19971 ; lAllison et al. 20"bg ). Allison 
et al. (2009) pointed out that the observed high-mass segre- 
gation can be of dynamical origin which can happen on very 
short time-scales of less than 1 Myr. The different degrees 
of segregation can be explained by the higher density of the 
ONC (about 15 times that of Tr 14) reducing the time-scales 



on which mass segregation occurs. In conclusion, a shallow 
PDMF as we observe in Tr 14 is not an uncommon feature 
but observed in several massive young clusters, but its origin 
remains elusive. 



5.3.2 'Traditional' Mass Function 

For comparison, we construct the PDMF also in the 'tra- 
ditional' way, counting stars in logarithmic mass intervals 
with Alog(m) = 0.2 dex and Alog(m) — 0.1 dex. Stars 
are counted in each bin and normalised to the considered 
area. The results are corrected for incompleteness and the 
resulting PDMFs are shown in Fig. 1101 Above 0.5 Mq, we 
fitted single-power law with the data points weighted by the 
individual Poisson errors. The corresponding slopes of are 
summarised in Table [2] The derived power law slopes agree 
within the errors with the PDMF of Fig. [9] although the 
PDMFs with 0.1 dex spacing appear steeper. Towards lower 
masses we also identify the drop of the mass function below 
0.5 Mq followed by a steep increase below 0.25 Mq. Due 
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Table 2. Slopes of the different mass functions 



Mass range Mq Binning r PD , P MS o"r PDiPM s r PD, a ii cr rpD aJ 

m c -3.2 1 star -0.46 0.20 -0.50 0.11 

0.1 dex -0.70 0.31 -0.60 0.26 

0.2 dex -0.59 0.27 -0.52 0.30 

0.25 - m c 1 star 0.73 0.59 0.63 0.32 
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Figure 10. MFs of Tr 14 derived with the 'traditional' 
method. Upper Panel: The PDMF through star counts 
with bins of Alog(m) = 0.2 dex. Above 0.5 Mq we 
find r PDiPMS = -0.59 ± 0.27 and r PDjall = -0.52 ± 0.30. Lower 
Panel: The PDMF obtained via star counts with a bin size 
of Alog(m) = 0.1 dex. Above 0.5 Mq we derive similar 
slopes as for the 0.2 dex-bins (r P D, P MSi = — 0-70 ± 0.31 and 
recall! = -0.60 ±0.26). 



to the sparsity of data points, we have not fitted the PDMF 
below 0.5 Mq. 



6 SUMMARY 

We presented deep H and Ks observations of the young 
massive cluster Tr 14, observed with MAD, the first MCAO 
system at the VLT. 

i) MAD (MCAO) in crowded regions: 
We revealed an impressive performance of the AO system 
over a large FoV. As part of the technical analysis we de- 
rived maps of Strehl ratios which show a low variation over 
the 68 arcsecx68 arcsec FoV. Mean Strehl values are mea- 
sured to 6.0 and 5.9 per cent in H- and 7\s-band, respec- 
tively. Simulating SCAO observations over the same FoV 
revealed smaller mean Strehls (1.6 per cent in H and 2.0 
per cent in Ks)- Although with lower maximum Strehl ra- 
tios (9.8 per cent in .ff-band, 12.6 per cent in JCs-band), 
when compared to SCAO systems, this shows a significant 
improvement of the spatial stability of the AO performance 



and demonstrates the opportunities offered by MCAO sys- 
tems for wide-field AO corrected observations. 

PSF photometry has been applied to derive prop- 
erties of Tr 14. The results show that the photometry 
is well es timated by usin g a second order variable PSF 
(see also ISana et al. 20101 ) . Photometric calibration was 
accomplished using seeing-limited NTT-SofT observations 
(|Ascenso et al. 2 007). A spatially constant ZP offset was 
found when compared to seeing-limited observations and 
used for photometric calibration. The resulting photomet- 
ric catalogue comprises 972 (1347) sources down to very- 
low mass (< 0.1 Mq) PMS stars. This doubles (triples) the 
number of detections when compared to seeing-limited ob- 
servations of the same field. 

ii) K s vs. H -K s CMD of TrU: 
The photometric catalogue with almost 1,000 sources was 
used to create a Ks vs. H — Ks CMD of Tr 14. It reveals 
a clear sequence of PMS stars, dozens of 'excess' sources 
and a sparsely populated MS. Comparing the CMD with 
PMS tracks by Siess et al. (2000), we find the CMD is best 
represented by an isochrone of 1 ± 0.5 Myr in age while 
hints for a somewhat older PMS population are observable. 
This places Tr 14 at a distance of 2.3 ± 0.4 kpc with a fore- 
ground extinction of Ak s = 0.38 ±0.03 mag (Av = 3.0 ±0.4 
mag). Hints for star formation over the last 3 Myr sup- 
port earlier claims of continuous star formation in the re- 
gion and a slightly older halo population of the cluster (e.g. 
lAscenso et al 20071 ). 

Hi) The LF and MF: 

The LF of Tr 14 is derived based on the Ks magnitudes. 
The observed turnover at Ks ~ 14.75 mag is consistent with 
a cluster of the derived age of ~ 1 Myr and the limiting mass 
at which Deuterium burning ceases. 

We derived the PMS stellar PDMF that appears to be 
well represented by a broken power law between 0.25 and 
3.2 Mq. We found a PDMF that shows a change of the 
power law slope at m c = 0.53 Mq comparable to that 
of a Kroupa-IMF. Such variations of the MF (flattening, 
turnover) of stellar systems are observed in other young 
clusters at similar stellar mass. The corresponding power 
law slopes, however, show a shallower PDMF when com- 
pared to the Kroupa-IMF (Salpeter slope F = —1.35). In 
the intermediate- and low-mass regime we find slopes of 
Ti = -0.50 ± 0.11 above m c and F 2 = 0.63 ± 0.32 below 
m c down to ~ 0.25 Mq. This reveals a deficiency in low- 
mass stars in the centre of Tr 14. Furthermore, comparison 
with other young massive clusters shows similarly shallow 
PDMFs in the centres of young and massive clusters. 

This study reveals the great improvement of AO correc- 
tions provided by MCAO systems in the case of wide field 
AO corrected observations. The stability of the correction 
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for atmospheric turbulence in combination with the correc- 
tion over a few arcminutes provides an ideal combination for 
investigations that depend on high spatial resolution instru- 
ments as well as a large FoV. Therefore, VLT-MAD gives a 
very promising impression of the future prospects of obser- 
vations with upcoming new telescopes and instruments. 
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